Spins coupled to a Z2-Regge lattice in 4d 
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We study an Ising spin system coupled to a fluctuating four-dimensional Z2-Regge lattice and compare with 
the results of the four-dimensional Ising model on a regular lattice. Particular emphasis is placed on the phase 
transition of the spin system and the associated critical exponents. We present results from finite-size scaling 
analyses of extensive Monte Carlo simulations which are consistent with mean-field predictions. 



1. INTRODUCTION 

Spin systems coupled to fluctuating manifolds 
are studied as a simple example for matter fields 
coupled to Euclidean quantum gravity. To de- 
scribe the gravity sector we used the Discrete 
Regge Model which is both structurally and 
computationally much simpler than the Stan- 
dard Regge Calculus with continuous link lengths. 
Here numerical simulations can be done more effl- 
ciently by implementing look-up tables and using 
the heat-bath algorithm. In the actual computa- 
tions we took the squared link lengths as = 
qi = 5/(1 + eai) with a; = ±1 and e = 0.0875. 
Because a four-dimensional Regge skeleton with 
equilateral simplices cannot be embedded in flat 
space, bi takes different values depending on the 
type of the edge I. In particular bi — 1, 2, 3, 4 for 
edges, face diagonals, body diagonals, and the hy- 
perbody diagonal of a hypercube. 

2. MODEL AND OBSERVABLES 

We investigated the partition function 
Z = J2 [ D[q]c^p[-I{q)-KE{q,s)], (1) 

where I[q) is the gravitational action, 
t i 
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The first sum runs over all products of trian- 
gle areas At times corresponding deficit angles 5t 
weighted by the gravitational coupling (5g. The 
second sum extends over the volumes Vi of the 4- 
simplices of the lattice and allows together with 
the cosmological constant A to set an overall scale 
in the action. The energy of Ising spins Si G Z2 , 

= (3) 

is defined as in two dimensions with the 
barycentric area Aij associated with a link lij, 
Aij = ^ , ^At- We chose the simple uniform 
measure as in the pure gravity simulations [Q, 
D[q] = Y[i dqiJ-{qi). The function ensures that 
only Euclidean link configurations are taken into 
account. 

For every Monte Carlo simulation run we 
recorded the time series of the energy density 
e = E/Nq and the magnetization density m — 
^jSi/Ao, with the lattice size Nq = L^. To ob- 
tain results for the various observables O at val- 
ues of the spin coupling K in an interval around 
the simulation point ii'o, we applied standard 
reweighting techniques 

With the help of the time series we compute 
the specific heat, C{K) = K^Noi{e^) - (ef), the 
(finite lattice) susceptibility, xi^) = No{{m'^) — 
{\m\y), the Binder parameter, Ul{K) = 1 — 
(m"')/3(r?T,^)^, and various derivatives of the 
magnetization, d{\m\)/dK, d\.T\{\m\) / dK , and 
d\n{m^)/dK. All these quantities exhibit in the 
infinite- volume limit singularities at Kc which are 
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shifted and rounded in finite systems. 

3. SIMULATION RESULTS 

In four dimensions it is generally accepted that 
the critical properties of the Ising model on a 
static lattice are given by mean-field theory, with 
logarithmic corrections. The finite-size formulas 
can be written as 101 



^ cx L(logL)s 
X (X (i(logi)I)^/^ 
K,{L) cx (i(logi)Ti)-i/^ 



(4) 
(5) 
(6) 



where the critical exponents of mean-field theory 
are a = 0, /3 = 1/2, 7 = 1, and i/ = 1/2. 

The gravitational degrees of freedom of the par- 
tition function (|l|) were updated with the heat- 
bath algorithm. For the Ising spins we employed 
the single-cluster algorithm [||. Between mea- 
surements we performed n = 10 Monte Carlo 
steps consisting of one lattice sweep to update the 
squared link lengths qij followed by two single- 
cluster flips to update the spins Si. 

The simulations were done for cosmological 
constant A = and gravitational coupling (3g = 
—4.665. This /3g-value corresponds to a phase 
transition of the pure Discrete Regge Model ||]. 
The lattice topology is given by triangulated tori 
of size Nq = L"* with L = 3 up to 10. From short 
test runs we estimated the location of the phase 
transition of the spin model and set the spin cou- 
pling Ko = 0.024 « Kc in the long runs. 

After an initial equilibration time we took 
about 100 000 measurements for each lattice size. 
Analyzing the time series we found integrated au- 
tocorrelation times for the energy and the magne- 
tization in the range of unity for all lattice sizes. 
The statistical errors were obtained by the stan- 
dard Jack-knife method using 50 blocks. 

Applying the reweighting technique we first 
determined the maxima of C, X; d{\m\)/dK, 
d\n{\m\) / dK , and d\n{m'^)/dK. The loca- 
tions of the maxima provide us with five se- 
quences of pseudo-transition points ifmax(i) for 
which the scaling variable x = (-ft'max(i) — 
K(,){L{\ogL)T^)~ should be constant. Using this 
fact we then have several possibilities to extract 



the critical exponent v from (linear) least-square 
fits of the FSS ansatz with logarithmic correc- 
tions (H), 

dUL/dK - {L{logL)^)^/''fo{x), (7) 
d\n{\m\P)/dK - (L(logL)A)i/-/p(x), (8) 

to the data at the various i^max(i) sequences. 
We also performed fits of a naive power-law FSS 
ansatz. The exponents l/v resulting from fits us- 
ing the data L = 4 — 10 are collected in Table |l|. 
Q denotes the standard goodness-of-fit parame- 
ter. For our simulations all exponent estimates 
with the logarithmic corrections and consequently 
also their weighted average l/i/ — 2.028(7) are 
in agreement with the mean-field value l/v = 2. 
With the naive power-law ansatz one also gets an 
estimate for l/iy close to the mean-field value, but 
clearly separated from it. 

Assuming therefore i/ = 0.5 we can obtain es- 
timates for Kc from linear least-square fits to the 
scaling behavior of the various -fCmax sequences, 
as shown in Fig. |l|. Using the fits with £ > 4, the 
combined estimate from the five sequences leads 
to Kc = 0.02464(4). 

Knowing the critical coupling we may recon- 
firm our estimates of l/i/ by evaluating the above 
quantities at Kc- As can be inspected in Table 0, 
the statistical errors of the FSS fits at Kc are sim- 
ilar to those using the iCmax sequences. However, 
here we have to take into account the uncertainty 
in our estimate of Kc- This error is computed by 
repeating the fits at Kc ± AiiTc and indicated in 

Table 1 

Fit results for l/f in the range L = 4 — 10 with 
a power-law ansatz with logarithmic corrections. 



fit type 




Q 


dU/dK at iiT^ax 


1.980(17) 


0.70 


dln{\m\)/dK at ifj^f^'"'^ 


2.032(10) 


0.59 


d\n{m^)/dK kI'^^""'^ 
weighted average 


2.038(10) 
2.028(7) 


0.55 


dU/dK at Kc 
dlii{\m\)/dK at Kc 
d\n{m^)/dK at Kc 
weighted average 


1.981(17)[13] 
2.027(9)[2] 
2.034(9) [2] 
2.025(6) 


0.70 
0.95 
0.85 


overall average 


2.026(5) 
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Figure 1. FSS extrapolations of pseudo-transition 
points -fCmax vs. (L(log L) 12 ) "^/"^ , assuming v = 
0.5. The error-weighted average of extrapolations 
to infinite size yields = 0.02464(4). 

Table |l| by the numbers in square brackets. In 
the computation of the weighted average we as- 
sume the two types of errors to be independent. 
As a result of this combined analysis we obtain 
strong evidence that the exponent v agrees with 
the mean- field value of = 1/2. 

To extract the critical exponent ratio ^jv we 
use the scaling (||) of the susceptibility x at 
its maximum as well as at i^c, yielding in the 
range L = 4 — 10 estimates of 7/1^ = 2.039(9) 
(Q = 0.42) and ^ijv = 2.036(7) [4] (Q = 0.85), 
respectively. These estimates for ^jv are consis- 
tent with the mean- field value of 7/t/ = 2. In 
Fig. ^ this is demonstrated graphically by com- 
paring the scaling of Xmax with a constrained one- 
parameter fit of the form Xmax — c(L(logL)4)^ 
with c = 4.006(10) (g = 0.17, L > 6). 

4. CONCLUSIONS 

We have performed a study of the Ising model 
coupled to fluctuating manifolds via Regge Cal- 
culus. Analyzing the Discrete Regge Model with 
two permissible edge lengths it turns out that the 
Ising transition shows the expected logarithmic 
corrections to the mean-field theory. We have 
also studied the pure Ising model on a rigid lat- 
tice without presenting the results in this short 
note. The critical exponents of the phase transi- 
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Figure 2. FSS of the susceptibility maxima Xmax- 
The exponent entering the curve is set to the 
mean- field value 7/1/ = 2 for regular static lat- 
tices. 

tion of the Ising spins on a static lattice as well as 
on a discrete Regge skeleton are both consistent 
with the exponents of mean- field theory, a = 0, 
P = 1/2, 7 = 1, and v = 1/2. In summary, from 
our comparative analysis with uniform computer 
codes we conclude that the phase transition of 
the Ising spin model coupled to a discrete Regge 
skeleton exhibits the same critical exponents and 
the same logarithmic corrections as on a static 
lattice. 
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